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ABSTRACT 



The formation of hydraulic jumps in a rotating, continuously 

stratified, Boussinesq fluid is investigated. A scale analysis indicates 

that jump-formation cases should be separated from non-jump cases 

2 

by the curve F = ARo where Ro and F are the initial Rossby and 
Froude numbers, respectively. Numerical solutions verify the curve 
up to F = 0. 3 and show that A is in the range 6.0 to 7,5. The solutions 
for F = 1, 0 are examined in considerable detail. The results of tliis 
paper extend the work of Houghton and Williams and Hori. 
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I. INTRODUCTION 



The formation of hydraulic jumps has been shown for a one-layer, 
homogeneous fluid for both non- rotating and rotating cases. This 
experiment will consider a one-layer, continuously stratified fluid by 
first investigating jumps in a non- rotating coordinate system and then 
extending the study to the rotating system. 

In a non- rotating one-layer fluid a hydraulic jump with a discon- 
tinuity in the velocity can arise from an initial state having any scale. 
If the system rotates, there may be a limiting initial scale beyond 
which jumps cannot form as shown by Houghton [1969] and Williams 
and Hcri [1970J. is '^xpnrrcci m ct continucu=>iv 

stratified atmosphere. In non-dimensional terms this corresponds 
to a critical Rossby number (Ro - U/fL, where U is the initial velocity 
scale, L is the initial space scale and f is the Coriolis parameter) 

below which jumps are excluded. This process also depends on the 

- 1 ^ 

Froude number (F = U/[D(g0Q d6o/dz) ], where D is the vertical scale 

and 6^ is the mean potential temperature). In the F-Ro plane there 

should exist a critical curve that would separate the jump cases from 

the non-jump cases as has been shown for a homogeneous fluid by 

Houghton [ 1969 ] and Williams and Hori [1970]. In this experiment it 

was shown that a critical curve exists and is valid for time periods of 
- ] ^ 

order H/[D(g0^ dO^/dz) ]. For longer time periods, however, jumps 



6 



m3.y occur in the non-jump region, indicating that the critical curve 
is actually a zone where the jump time increases rapidly. 

There are several possible applications of this work in meteorology 
and oceanography. Rotation does have an effect on mesoscale motions 
in the atmosphere, but the influence is not dominant as it is in large- 
scale motions. Tepper [1952] has proposed that squall lines are 
modified hydraulic jumps. Freeman [1948], Abdullah [I 966 ], and 
Houghton and Kasahara [I 968 ] have proposed other mesoscale hydraulic 
analogies. Sasaki [1959] and Ogura and Charney [1962] have studied 
squall lines in continuously stratified atmospheres. 

Internal waves in the ocean are affected by rotation, and the break- 
ing of th‘ ► o -vr p 3 3 ould be described by this model. Rapidly moving 

frontlike temperature disturbances have been observed in the ocean 
[Gowans, 1969] and these may be internal waves modified by rotation. 

In this experiment the Boussinesq equations are given and the 
modeling relations are introduced in II. In the model the time-dependent 
quantities are functions of x and z only. The quantity f is assumed to 
be constant for each independent case study. The finite difference 
approximations are presented in III. A scale analysis is presented 
in IV which is useful in the discussion of the numerical solutions. The 
initial conditions are given in V and tlie numerical solutions are pres- 
ented and discussed in VI. The initial state consists of a sinusoidal 
disturbance (amplitude dependent on F) which is superimposed on the 
basic state. After a sufficient period of time, the sine wave is distorted 
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I 



in such a way that a hydraulic jump is formed. In VII conclusions and 



recommendations for further investigation are given. 



8 



II. BASIC EQUATIONS 



In this experiment the Boussinesq equations are employed and 
the domain is bounded by two rigid horizontal planes. The compress- 
ibility of the atmosphere is neglected in the Bossinesq approximation. 
This should not be of qualitative importance since the density scale 
height of the atmosphere is much larger than the vertical scale of the 
typical jump zone. The hydrostatic approximation is also used in 
this experiment and it should be good as long as the horizontal scale 
is much larger than the vertical scale [Lamb, 1945, p. 258]. 

The Boussinesq equations which were obtained by Ogura and 
i’iiillipci vn tlie 

rotation is included: 




(II. 1) 



If + v-ev 




(II. 2) 



'7-V + IT = 0 and 



(II. 3) 



^0 3 ® 






(II. 4) 
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where: 



0Q = constant reference potential temperature, 

P = constant reference pressure, 

« = R/Cp. 

6 = T(Pq/P) - ~ departure of the potential temperature 

from 0Q, and 

0 = Cp^o'P/Po)" + S^- 

Heating and friction have been neglected in these equations even though 
they may be important under certain conditions. 

All quantities are independent of y in this model. The Coriolis 
parameter, f, is assumed to be constant throughout the model but is 
varied for individual cases. A Cartesian coordinate system is intro- 
duced with X increasing toward the east and y increasing toward the 
north. The x and y velocity components are u and v respectively. 

For this model, Equations (II. 1)-(II. 3) can be written in the 
following form: 

. ^(uu) , i(uw) , r n (II. 5) 

Tt* + TV TV" T7 - - 0, 

AV ^ A(WV) ^ A(VW) X., _ n nr 

H + -nr ^ rr + 







jf«e) 


(II. 7) 


At 


1 


Ax AZ 




A u 




Aw ^ 


(II. 8) 


Ix 


*!• 


IT " 0. 
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The boundary conditions are 



w = 0, z = 0 , D (II. 9) 

where D is the distance between the horizontal plates. 

The integral of Equation (II. 4) gives the departure of 0 from its 
vertical average, and this average can be obtained as follows. The 
following averages are defined: 

O 

where 2 n L is the periodicity of the fields in x. If Equation (II. 8) is 
averaged vertically and (II. 9) employed 

»- <«>=<(“)>• 

i X 

Similarly, the vertical average of (II. 5) is 






(II. 12) 



If this equation is in turn averaged with x and if periodicity in x is 
imposed, Equation (II. 12) can be written 



d<(u)> 



cJi 



- = 0 . 



(II. 13) 



The quantity ^(V)^ represents the net mass flux in the y direction. If 
<(V1> were different from zero, ((wP could not be zero if a wall 
were present at some value of y. Since this would violate the boundary 
conditions [Equation (II. 9)], the following is a requirement in this 
experiment: 
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<(v)> = 0. 



(II. 14) 



It then follows directly from (II. 13) that 

d<(u)> 



d't 



0 . 



(II. 15) 



Now from (II. II) 
reduces (II. 12) to 



><u> d<(w)> 









which together with (II. 15) 

+ (II. 16) 



The corresponding equation is obtained by averaging Equation (II. 6) in 
X and z and by introducing (II. 14), i. e. , 

Equation (II. 17) shows that the x- and z-averaged zonal wind is zero. 

If E4uaUons (II. I), ;II. 16} - -id (IL 17} nzra, Enua^-rcns {II. S) 

(II. 6) can be rewritten in the forms 



jyK 


•f 


> ( uw) 


^ i(uw) 








>1 










it 




i ( w v) 


, i(wv) r 



(II. 18) 



(II. 19) 



By integrating the hydrostatic equation [Equation (II. 4)] the following 
expression for is obtained: 








(II. 20) 



The following energy integral is an invariant for this model: 



£• = [<( «')) * (( 




( 11 . 21 ) 
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c. 



III. FINITE DIFFERENCE APPROXIMATIONS 



Equations (II. 7), (II. 8), (II. 18) and (II. 19) were solved numerically 
by introducing finite differences in x, z and t. The solutions were 
obtained by a simple marching process, but if the horizontal variations 
were two-dimensional, Poisson's equation would have to be solved at 
each time step [Smagorinski, 1958; Hinkelmann, 1959]. 

The vertical domain was divided into 2K layers of equal thickness, 
and the surfaces separating the layers were denoted by k. The arrange- 
ment of variables on the bounding surfaces is shown in Figure 1. The 
divergence of the vertical flux of the scalar S at z = Zj^ is approximated 

cf s 



i(wS)' 



k 



(S^.i ~ Wfc.' (St 

2 -Zk-P 



(III. 1) 



,, K JLJ±LL1J..1 UdJJJJJUJlL I . D 

k.K-i 





K - 2 

kH 

k.i 

k*0 



yv 

yz^J 

w 

Trrnmhiiiijiiii 



Fig. 1. Arrangement of variables in z. 
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The finite difference approximation to the integral in Equation (11.20) 
used is 



=Z (V.-Zj)] +C , l<>2 • 



rZl 

The value of C = is immaterial since it will disappear when the 

vertical average is removed. The finite difference integral of Equation 



(II. 8) used is 






(III. 3) 



The vertical grid structure was introduced by Lorenz [i960] for use in 
the balance equations. 

All of the dependent variables occur at each x grid point. The x 
derivatives were approximated as 



/ il) jj ~ S;. 
Wx/i “ 2Ax 



3nd 



(III. 4) 



(u.„ -t- U/)(S, + I +S;) - (u,- + + 



[^ 1 - 



4 Ax 



(III. 5] 



Here i denotes the grid index and Ax the mesh length. This is a variant 
of the finite difference scheme that Arakawa [I966] developed [see also 
Lilly, 1965 ]. The Arakawa scheme is designed to avoid the non-linear 
computational instability which was discovered by Phillips [1959]. In 
particular, the finite difference approximation to the integral, 
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i 






(III. 6) 



is zero when (III. 1) and (III. 5) are used. The integral I is always 
zero in the continuous case. If the domain averages of u , v and 6^ 
are conserved by the non-linear terms throughout a numerical integra- 
tion, then no non-linear computational instability can occur. When 
finite differences in time are employed, these quantities will not be 
identically conserved and instability may arise. 

The continuous equations conserve the total energy E, which is 
given by Equation (11.21) , If non-linear computational instability 
occurs, E cannot be conserved. It was found that E remained nearly 
constant in the numerical solutions in this experiment, which indicates 
that the spurious computational production of energy was small. This 
was important because it was desired that the computational scheme 
would not introduce a small-scale energy source which could produce 
a false jump. 

Centered time differences were used throughout this experiment 
with the exception of the first time step which was a forward time step. 
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IV. SCALE ANALYSIS 



A scale analysis is now performed that will indicate the dependence 
of the solutions on the Rossby number and Froude number. The inde- 
pendent variables are written 

t = Tt', X = Lx' and z = Dz'. (IV. 1) 



The dependent variables are scaled as follows: 

u = Un v = Vv’, w=Ww‘, 0 = UDN0 ' and 
6 = (UN0o/g)9’ . (IV. 2) 

The Rossby number and Froude number are defined as 

Ro = U/fL and F = U/DN (IV. 3) 




The time scale T and the v-scale V are functions of Ro and F while 
the w-scale is W = UD/L. 

An examination of the case where F;^ 1 and F^ Ro is now given. 
This should parallel the results obtained by Houghton [1969]. The 
appropriate time scale is T = L/DN, and the v-scale is V = (F/Ro)U. 
The later was obtained following Houghton by balancing with fu 

in Equation (II. 6), The non-dimensional equations for this case are 



iV ^ hy' ^Z' 





(IV. 5) 
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(IV. 6) 



^v' 


t f 


Mu'v') 


i(v'w') 


Ji' 




^2' _ 




+ F 


' ^(u'e') 


d (w'0 ) 


H' 


'hr 


^2' 



+ U - 0, 



+ w' = 0, 



(IV. 7) 



I? ■'■ T? *° ®'’‘' 



(IV. 8) 



ii' - 9' ^ 0 . 

H' 



(IV. 9) 



Hyciraulic jumps are formed through the action of the non-linear terms 
c\ /T‘\r \ --..ri My \ I'yvr.n when uic 

small (F« 1), they can produce a jump in a non-rotating system. If 
a hydraulic jump is to be prevented, the Coriolis term in (IV. 5) must 
dominate the non-linear advective terms, which implies that 



pi 



(IV. 10) 



The curve dividing the jump region from the no-jump region should 
then have the form 

where the constant A should be appreciably larger than 1. Houghton's 
analytical curve corresponds io AOi 6.5; Williams and Hori [1970] 
found that A = 6. 0 to 7. 5 was in agreement with their numerical 



17 



solutions. As in Williams* and Hori's experiment, with A in the range 
6. 0 to 7. 5, the above scale analysis does not apply for F Ay 1, since F 
is then greater than Ro. 

Two regions in Ro are of interest when F a^ 1. If Ro a^ 1, all 
terms in Equations (IV. 5), (IV. 6) and (IV. 7) are of the same order 
and a jump would be expected to form. However, when Ro 1, the 
appropriate time scale is T = f " ^ and the v- scale is V = U. The non- 
dimensional equations then become 




(IV. 12) 



(IV. 13) 



(IV. 14) 



along with Equations (IV. 8) and (IV. 9)- If the terms in Ro are neglected, 

the resulting equations describe an inertial oscillation in u and v; the 

temperature field is independent of time. 

The scale analysis indicates that, when F is small, the jump 

2 

boundary is given by F = ARo , and that as F approaches 1, the 
boundary will be between Ro =1.0 and Ro = 0. 1. These results were 
tested by direct numerical integration of the basic equations. 
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V. INITIAL CONDITIONS 



\0 

The initial x-averaged stability ~ is independent of z. A wave 
sinusoidal in x, the amplitude of which is directly dependent on the 
Froude number, is superimposed upon this basic state. 

The initial fields for this experiment are given by the equations 

w(^,0) = U cos(kx) |cos[-| (z -|)] ' :|j- , 



v(x,0)- 0 and 



(V. 1) 



codloc) 



3 

The constants which define the initial state have the following 
numerical values: 

D - 

■ (V.2) 



If = 4.0 ‘'K kv>r‘ and 



1 - 0.0526*5 m sec^ ®K‘\ 

The wavelength of the disturbance, L, is equal to 1, 000 km. 

After finite differences are introduced, it is necessary to adjust 
the values of u which are given by (V, 1) in order that 




(V.3) 



This merely states that the finite difference vertical velocity at the 
upper boundary must be zero. 
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The following initial conditions were also examined: 

u (x^o) = U cos(kx)cOS (^) , 

v(x,0)s0 5nd (V.4) 

6 (,, 0 }..^ cos(b)sin(f) + i(f-2). 

When f = 0 they lead to an exact finite amplitude wave solution 
which moves with the linear gravity wave speed. In this solution the 
non-linear terms in Equations (II. 7) and (II. 18) cancel out identically. 
In this case no jumps can occur although they were observed for f 5 ^ 0. 
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VI. NUMERICAL SOLUTIONS 



The numerical solutions are presented in non-dimensional form. 
The time for all results presented is given in units of the time scale, 
L/DN. Most of the results are for 1. Enough results are pre- 

sented for F 4^ 1 in order to determine the jump curve. This will 
provide a basis for comparison with the work done previotisly by 
Houghton [ 1969 ] and Williams and Hori [1970] in the one-layer, homo- 
geneous system. 

A quantity useful in the analysis of the jump criterion is defined 
as 




The quantity is approximated by (u^^j - u^)/ A x. Here the sub- 
scripts max and min refer to the maximum and minimum overall 
points at a given time. The non-dimexisional quantity d is a measure 
of the width of the zone through which u varies most rapidly. It will 
be an overestimate of this width when the values of u at the edges of 
the zone are not near u ^ and u . However, if a discontinuity 
such as a pressure jump does occur, d must approach zero, since it 
is inversely proportional to the gradient of u. 

Figure 2 shows the variation of d as a function of time (time unit 
L/DN) for F = 1. 0 and Ro = 00 . Results are shown for three values of 
A x = 2 7 T L/n corresponding to n = 2-0, 40 and 80. As can be seen, d 
decreases to a minimum at approximately DNt/L = 0.4. 
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Fig. 2. The quantity d as a funct on of DNt/L for F-1.0 and Ro- 
Each curve corresponds to a iifferent value of x=2 L/n. 
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Detailed solutions are presented for three cases in which F = 1. 

In these numerical solutions A x = 2 r L/n, where n = 80. Figure 3 
shows the time evolution of the non-dimensional variables u* = u/U 
v' = v/U and 0* = (0 - 0) /(UN0^/g) [where 9 represents the horizontal 
average] for Ro = 1. The u and 0 fields show that a hydraulic jump is 
forming at DNt/L = 0. 37. The v field, v/hich begins at zero, grows in 
opposite phase to the u and 0 fields through the action of the Coriolis 
force. Once the v field is present, the Coriolis force will tend to 
reduce the amplitvxde of the u field while the amplitude of the 0 field 
remains relatively unchanged throughout the period of integration. For 
this value of the Rossby number the Coriolis force was unable to 
prevent the formation of the hydraulic jump. 

Figure 4 contains the evolution of the non-dimcn&lonal variables 
u*, V* and 9* for Ro = 0. 1. It appears that an inertial oscillation is 
present with the 0 field remaining very nearly constant throughout the 
integration period. It was found that when the integration time was 
extended, though, the fields never were restored to the original as 
was shown for the one-layer, homogeneous case by Williams and 
Hori [1970]. 

The results for Ro = 0. 35 are presented in Figure 5 using the 
same format as above. In this case there is no evidence of any 
reduction in the u field as the jump time is approached. 
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Fig. 3. The fields u’, v’ and 
0' at vario\is times for F=l. 0 
and Ro = 1 . 0. 




Fig, 4. The fields u', and 
0' at various times for F=l. 0 
and Ro=0 . 1 . 
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1.5 




and Ro=0. 35 . 
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Results for F = 0. 04 and Ro = 0. 09 are contained in Figure 6 
showing the evolution of u', v‘ and 0’ for a typical jump case when 
F 1. Figure 7 contains the evolution of the variables for F = 0. 04 
and Ro = 0. 07 giving the results for a non-jump case when F is again 
much less than 1. 

In Figure 8 and Figure 9 the results are presented for F = 0. 3 
and Ro = 0. 3 and 0. 1, respectively, showing the time evolution of 
u‘, V* and 0' when F is in the range 0. 1 to 1. 0. 

The above figures are taken from the results of the lowest layer 
in the model. Figure 10 is included in order to show what is taking 
place in the other layers at time zero and at an approximate jump 
time. The case presented is for F = 1. 0 and Ro = 1. 0. As can be 
seen, only the lower half of the model is shown. This provides a 
picture of what is taking place in the most active of the layers in the 
model. 

The d procedure was used to determine the curve separating the 

jump cases from the non- jump cases over the range F = 0. 04 up to 

F = 1, 0. These results are shown in Figure 11, with the logarithm 

2 

of F plotted against the logarithm of Ro. The lines F = 6. ORo and 

2 

F = 7. 5Ro are included for convenience. These results for Equation 
{IV. 11), which was derived by scale analysis leading to the conclusion 
that A is in the range 6. 0 to 7.5. This relation does not hold for Fr^ 1, 
but this was to be expected since the analysis required that F ^ Ro. 
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J 



0 



I « 




and Ro = 0. 09 . and Ro = 0.07. 
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r 





o' at various times for F=0. 3 O' at various times for F=0. 3 

and Ro=0.3. and Ro=0. 1. 
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Fig, 10. The lowest five layers of the model at time zero and at 
approximate jump time for F=l, 0 and Ro=l. 0. 
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$ 




Fig. 11. Summary of jump-formation experiments with the logarithm of 
F plotted against the logarithm of Ro. The crosses indicate jump 
occurrence and the dots indicate cases where no jumps were found. 
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The non-jump solutions included in Figure 11 are determined over 
an integration time of L/DN. 
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VII, CONCLUSIONS AND RECOMMENDATIONS 



The formation of hydraulic jumps in a rotating, continuously 

stratified fluid has been investigated. A scale analysis indicated that 

jump cases should be separated from non~jump cases by the curve 
2 

F = ARo in the F-Ro plane. The Rossby number Ro and the Froude 
number F were utilized as direct inputs to the equations and used as 
a means of determining the initial scales (U = FDN). It was determined 
that the results of this experiment (with F appropriately redefined) 
closely paralleled the results of Houghton [1969] and Williams and 
Hori [ 1970 ] in which one-layer, homogeneous fluids were studied. 



range 6. 0'7. 5 and that the relation F = ARo^ holds up to at least 
F = 0. 3. Extrapolation to F = 1, 0 showed jumps occurring for 
smaller values of Ro than were given by the above curve. Thus it can 
be seen that the jump criterion becomes more independent of F as 
F = 1, 0 is approached. 

This work should be extended to an atmosphere in which the jump 
region is bounded by a stable layer rather than a rigid plane as in 
this case. Such a study could require additional layers with further 
consideration given to layers of varying thickness. Further considera. 
tion should be given to moisture in future experiments. 
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